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Abstract 



Tij" ■ Dispersion curves to a oscillatory reaction-diffusion system with the self- 

(N, 

consistent flow have obtained by means of numerical calculations. The flow 

V ' 

(•^ ■ modulates the shape of dispersion curves and characteristics of traveling 

^^ ■ waves. The point of inflection which separates the dispersion curves into 

o 
o\ 



two branches corresponding to trigger and phase waves, moves according to 
the value of the advection constant. The dynamics of phase wave in reaction- 
diffusion-advection equations has been studied by limit cycle perturbations. 



^ ■ The dispersion relation obtained from the phase equation shows that the com- 

petition between diffusion and advection constants modulates the oscillation 
frequency from the bulk oscillation in the long-wave dynamics. Such a compe- 
tition implies that phase waves with the flow have a wider variety of dynamics 

than waves without the flow. 
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I. INTRODUCTION 

In biological systems of cell population, taxis induces cooperative movement of cells 
or the flow of cellular mass. Cells behave according to chemical and physiological signals 
given by themselves. Two well-known pattern formations for population dynamics of cells 
are the aggregation of cellular slime molds and the branching growth of bacterial colonies. 
Slime mould cells secrete signaling chemicals responding to extracellular chemicals, and 
move towards increasing some chemical concentration. They aggregate to centers of the 
wave pattern of chemicals with dendritic streaming lines and eventually form slugs [|T],0,^. 
Morphology of bacterial colonies depends on agar concentration and nutrient level. Bacteria 
are swimming toward high concentration of nutrients, and form various colonial patterns as 
concentric rings and dendritic branchings @,§,|1. Mechanisms of these pattern formations 
have been modeled by taxis equations or reaction-diffusion systems with the cellular flow. 

For one of protozoan myxomycete, the Physarum plasmodium, which is a unicellular 
organism, a reaction-diffusion-flow model is presented in relation to control mechanism 
of amoeboid movements. The flow of endoplasm and chemicals induced by protoplasmic 
streaming in Physarum plays an important role in individual behavior under a large Plas- 
modium state. It is not difficult to look for reaction, diffusion and flow couplings in various 
biological process, more and more. The self-consistent flow as above is realized by biolog- 
ical active process, and is a different character from physico-chemical ones. Effects of the 
self-consistent flow on simple reaction-diffusion systems are to be clarified in order to study 
a mechanism of functional self-organization of multicellular systems. 

In the present article, we discuss basic effects of the flow on oscillatory phenomena in a 
reaction-diffusion system with the self-consistent flow. Although we dealt with a model for 
the Physarum plasmodium, it is expected that the results obtained here are applicable to 
other biological systems with the flow. We carry out numerical calculation of the dispersion 
relation in the system, and show effects of the flow on traveling plane waves. Then we derive 
the dynamics of phase waves from reaction-diffusion-advection equations in general form 



by means of limit cycle perturbations. In the coefficient of the gradients of the phase, the 
advection terms compete with diffusion terms. 

II. MODEL EQUATIONS FOR THE PLASMODIUM 

The Plasmodium of Physarum has a cytoplasmic cortex (ectoplasm) ffiled with endo- 
plasm. The ectoplasm shows contractile oscillation everywhere within the organism [0,^,^, 



and the contraction causes intracellular streaming of the endoplasm ||I^. The contraction- 



relaxation behavior is regulated by metabolic cycles of chemical oscillators [0,0. Metabolic 
chemicals transported by endoplasmic streaming, affect chemical oscillators ||T3| , p!^ , p!5| . 

A model of contractile and motile dynamics of the Physarum plasmodium is represented 
by conservation of the cytoplasmic mass and reaction-diffusion system of metabolic elements 
with endoplasmic streaming [|1^,0]. Under some assumptions, oscillatory dynamics in the 



Plasmodium is governed by the following reaction-diffusion equations with advection terms: 

-^ + MVu ■ Vu = f (u) + Z^V^u, (1) 

at 

where u is a A^-component vector of reaction species which are separated into two types 

of metabolic elements, free chemicals and non-free ones. While the former is transported 

by the protoplasmic streaming, the latter is bound or stored at some cytoplasmic structure. 

Reaction kinetics f denotes the metabolic oscillation which emerges from the coupling of two 

types of reaction elements mentioned above, and the system is assumed to have a limit cycle 

orbit. The quantity Z) is a diagonal matrix of diffusion constants. A tensor M represents 

advection coefficients, here MVu is a flow vector for each reaction component. This system 

has a generic form of react ion- diffusion equations with the self-consistent flow, and we derive 

the phase dynamics from eq. |I] in section |V. 



In our numerical calculations, we use a more simple two-variable model [|T7 

^ = ^K^), (2) 



where u and v are the concentrations of a free chemical substance and a bound/stored one, 
respectively. We adopt the Schnackenberg's tri-molecular two species model for reaction 
kinetics of chemicals [|I^, 

f{u,v) = a — u + u'^v, g{u,v) = b — u^v, (3) 

here a and b are positive constants. In the spatially homogeneous conditions {D = and 
w = 0), the system has a stable limit cycle as shown in Fig. |l| for b — a > (a + b)^. The 
quantity D is the diffusion constant of the free chemical u. The velocity of the endoplasmic 
flow w is determined by the concentration of the metabolic chemical as 

w = qVu, (4) 

and thus the flow is self-consistent. Here q depends on the permeability and the mechanism 
of intracellular pressure. In the following calculation, q is assumed to be constant. 

III. DISPERSION RELATION 

The numerical calculations were carried out for eq. ^ on a ring, that is a one-dimensional 
region with periodic boundary conditions. Using the limit cycle, we initiate a pulse traveling 
on the ring, and solve eq. ^ until the solution becomes periodic in time. After we measured 
the rotating period of the traveling pulse on the ring, we repeat the calculation for rings 
of different lengths. As above, the dispersion relation has been obtained for periodic wave 
trains with the stable propagation. 

In the measurements, parameters have been set as a = 0.1, b = 0.5, and D = 1.0. We 
have used the explicit Euler's method for reaction terms, upwind differencing method for 
advection terms, and the implicit method for diffusion terms. The data have been recorded 
for waves in stationary propagation. 

Dispersion curves in the reaction-diffusion-advection system ^ are shown in Fig. ^ (wave 
number vs frequency) and Fig. |^ (period vs velocity) for various values of the advection 



constant q. Variations in the self- consistent flow make the propagation of plane waves quite 
changed. The inflection point which separates two branches corresponding to trigger waves 
and phase waves |l9| , p0| , moves according to the flow. It is the point that the variation of 
propagating behavior depends on the wave length. The other point is that the sign of q 
gives different effects on waves in propagation features. 

For the positive value of g, the oscillation frequency uo of phase waves gets slightly greater, 
while that of trigger waves becomes smaller. Hence, phase waves travel a longer way for the 
greater value of q even if they have the same frequency. Such a case of positive q corresponds 
to the phenomena of phase waves in the Plasmodium, when we regard the reaction kinetics of 
eq. ^ as the oscillation of some metabolic element like Ca^+. The flow make the Plasmodium 
communicate local information to the wide area in the cell with phase waves. 

When the advection constant q is negative, the opposite situation to the case of positive 
q arises. Furthermore, we find out the steep rise of the frequency in Fig. ||. Such steepness 
denotes the sharp transition from phase waves to trigger waves with increasing dimensions 
of the propagating media. 

IV. PHASE DYNAMICS 

In this section, we show that the effect of the flow comes out through the coefficient of 
the nonlinear term in phase dynamics for reaction-diffusion-advection systems. By means 
of limit cycle perturbations, the dynamics of phase waves in ordinary reaction-diffusion 
systems for oscillatory media are described by the Burgers equation pT] , p2| . We adopt the 
similar method to oscillatory reaction-diffusion equations with advection terms introduced 
in section |I[ 

We assume that the limit cycle of eq. |l| has the frequency uoq. Then a solution of homo- 
geneous oscillation is 

u = Uo(r), r = uj^t, where 

ujqVLq = f(uo), and Uo(r + 2tc) = Uo(r). 



Since the system |T| is invariant under the time translation, Uo(r + ip) {ip is an arbitrary 
constant) is also a solution of eq. |I]. 

Let us consider the dynamics of phase waves in eq. |1| when the advection and diffusion 
terms are small. We introduce multiple scales and asymptotic expansions, 

R = yer, r = uqI, T = et, 

u = Uo(r + ^) + eu H , (5) 

where e is a small parameter and ip = i^iR^T). Substitution of eq. |^ into eq. |l| yields 
perturbation equations for each order in e: 

d di 

Luj = hj, L = uJo—-—{uo), (6) 

here j = 1,2, . . ., and bj denotes the inhomogeneous term of the jth order equation. For 
the first order equation in eq. ^ the inhomogeneous term is 

bi = -^ - MVrUo ■ VrUo + DVluo 
= -u'o^-MuX|Vi,^r 

where Vr is the nabla operator in respect to scaled coordinates R. Thus the solvability 
condition for ui gives the dynamics of phase waves: 

|^ = ciV^V' + C2|V^^p. (7) 

If we use a new dependent variable V = VriP, then it satisfies the Burgers equation. The 
coefficients Ci and C2 are obtained from the relations, 



Cj = (v^vJ)/(v^uo), 

vi = Du'o, V2 = Du'^ - Mu[,u: 



0' 



here (v"!", v) = / d^(v^, v) and v'l" is the nontrivial periodic solution to the adjoint dif- 
ferential equation of Lv = 0. Equation |^ describes slow and slight modulation of the 
homogeneous oscillation with the frequency ojq by the phase ■?/'. We note that the coefficient 
of the nonlinear term, C2, show competition between diffusion and advection. 
When we use the quantity = u^t + ■?/', eq. becomes 

^ = ^o + ciV'0 + C2|V0|^ (8) 

The dispersion relation is thus estimated from the phase equation ^ through the wave char- 



acteristics uj = d(f)/dt and k = V0 as ||T9 



u = Uq + C2k'^ + ■ ■ ■ , k = \k\. (9) 

Since the scaling of coordinates in the perturbation expansions ^ means spatially slight 
modulation, k = 0{^/e), eq. ^ is the Taylor expansion for the dispersion curve, uj = uj{k) 
in the vicinity of A; = 0. Thus, the coefficient of nonlinear term in eq. ^ is C2 = u;"(0)/2. 
Here eq. |^ has no linear term in k because of the reflectional symmetry in space of eq. |1|. 
As mentioned above, C2 depends on advection constants as well as diffusion constants, and 
hence C2 can take the value of a wide range. 

We point out that the dispersion relation is only applicable to some of periodic waves 



with stationary traveling. It is not adopted to waves with non-uniform phase gradients pO 
In such a case, we need to use the phase equation ||, or analyze eq. |I| directly. 

V. DISCUSSION 



By means of numerical calculations for eq. ^, it has been shown that the flow make 
remarkable differences in the dispersion curves or propagation of phase waves. To elucidate 
such differences induced by the flow, we consider the spatial scale. The two-component 
system ^ has the diffusion term only in the free element. Thus, without advection term, 
the variation of the diffusion constant is canceled by the scaling of spatial coordinates. In 



contrast, such a cancellation is impossible with the advection term because the advection 
coefficient is independent of the diffusion constant. 

The deformation of the dispersion curves with the flow term has been shown for the 
branch of phase waves by the phase equation derived from eq. 0. Furthermore, we have 
found out another phenomenon of phase wave mentioned below. 

Advection coefficients govern the coefficient C2 of the quadratic term in the wave number 
k to the dispersion relation ^ The coefficient C2 is given by the addition of two parts which 
stem from the diffusion and advection terms, respectively. Thus, the competition between 
the diffusion and advection constants governs the value of C2- The variation of C2 gives a 
quadratic change in the frequency of long waves. This result elucidates the branch of phase 



waves {k ~ 0) for the dispersion curves obtained numerically in section [III. 

Varying the advection constants, we can take C2 negative. Such a situation seems to 
be impossible for simple reaction-diffusion equations without the advection terms. The 
dispersion curve is convex for C2 < 0, and hence cj = at some wave number. This means 
that the phase wave is frozen and stationary. It sounds strange, and we need to study these 
waves in detail. One of such a system is the X-u system with advection terms of phase 
gradients [P5|. 



When the coefficient C2 becomes sufficiently small, it may be possible to derive other types 
of phase equations by introducing different scaling of independent variables from ones used 
in the present analysis. Furthermore we need to study waves with relaxation oscillations 
to know the characteristics of trigger waves with the flow, because the phase equation is 
satisfied only in the region of the phase-wave branch. 
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FIG. 1. The limit cycle of the Schnackenberg's tri-molecular two species model. 
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FIG. 2. Dispersion curves obtained by numerical calculations. The oscillation frequency cv 
monotonically increase with the wave number k, and has the bulk frequency ujq = 0.495 in the 
limit of /c ^ 0. 
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FIG. 3. The propagating velocity v is determined as tlie pliase velocity to/k. The velocity 
becomes large and goes to infinity at the bulk period tq = 12.7. For negative value of q, the curves 
have large convex regions in period corresponding to trigger waves. 
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